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Abstract 


Continental shelves are believed to play a major role in carbon cycling due to their high productivity. 
Particulate organic carbon (POC) burial has been included in models as a carbon sink, but we show 
here that seasonally produced dissolved organic carbon (DOC) on the shelf can be exported to the 
open ocean by horizontal transport at similar rates (1-2 mol C m‘ yr ) in the southern U.S. Mid- 
Atlantic Bight (MAB). The dissolved organic matter (DOM) model imbedded in a coupled 
circulation-biogeochemical model reveals a double dynamics: the progressive release of dissolved 
organic nitrogen (DON) in the upper layer during summer increases the regenerated primary 
production by 30 to 300%, which, in turns, enhances the DOC production mainly from 
phytoplankton exudation in the upper layer and solubilization of particulate organic matter (POM) 
deeper in the water column. This analysis suggests that DOM is a key element for better representing 
the ecosystem functioning and organic fluxes in models because DOM (1) is a major organic pool 
directly related to primary production, (2) decouples partially the carbon and nitrogen cycles 
(through carbon excess uptake, POM solubilization and DOM mineralization) and (3) is intimately 
linked to the residence time of water masses for its distribution and export. 
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1. Introduction 


Dissolved organic matter (DOM) has received increasing attention over the past few decades 
because dissolved organic carbon (DOC) represents by far the largest pool of organic carbon in the 
ocean. DOC export from the surface global ocean is estimated at 20% of total organic carbon flux to 
the deep ocean (Six and Maier-Reimer, 1996; Hansell, 2002), which represents an important control 
on atmospheric carbon dioxide levels (i.e., the biological pump). Particulate organic carbon (POC), 
which accounts for 80% of the organic carbon export in the open ocean, is mainly recycled on the 
shelf and fuels the DOC pool. If residence times of shelf waters are shorter than the lifetime of much 
of the seasonally produced DOC, the horizontal DOC flux could represent the main contribution to 
the export of organic carbon to the open ocean and exceed sinking POC fluxes (Bauer and Druffel, 
1998). 

The Mid- Atlantic Bight (MAB) is the central region of the eastern U.S. continental shelf 
characterized by high rates of primary productivity and strong residual circulation. This region is 
thus a potential area of organic carbon export to the open ocean. The DOC pool in the MAB is one 
to three orders of magnitude larger than the POC pool (Bauer, 2001). The hydrography and 
circulation of the MAB is well studied (Biscaye et al., 1994) with a general north to south flow from 
Georges Bank to Cape Hatteras. The greatest inflow of water to the MAB (0.4 Sv, Beardsley and 
Boicourt, 1981) is from Georges Bank and is characterized by relatively low temperature and 
salinity. This flow represents the southern extension of the Labrador Current with averaged (total) 
DOC concentrations between 72 and 75 mmol m" (Vlahos et al., 2002). In the southern MAB near 
Cape Hatteras, the warm and salty water of the Gulf Stream has a major impact on the flow and 
exchange of water on the shelf and slope. Although the mean circulation is along-shelf in the 
southwestward direction, cross-shelf exchanges due to frontal instabilities and eddies displace 
significant portions of the shelf water to the open ocean north of Cape Hatteras (Biscaye et al., 
1988). The freshwater inflow (0.005 Sv) to the MAB is about 1% of the total water inflow 
(Beardsley and Boicourt, 1981), but represents a significant input of total DOC with mean 
concentrations of 200 to 400 mmol C m" in the mid-Bay of the Chesapeake Bay (data from the 
Chesapeake Bay Program, http://www.chesapeakebay.net/wquality.htm) throughout the year and 90 
to 190 mmol Cm' 3 in the mouths of the MAB estuaries (Bates and Hansell, 1999; Fisher, 1998; 
Harvey and Mannino, 2001; Sharp et al., 1982; Taylor, 2003). Terrestrial DOC fluxes and 
composition in the MAB region were investigated (Delaware Bay, Mannino and Harvey, 2000a; b; 
Chesapeake Bay, Mitra et al., 2000; MAB, Aluwihare et al., 2002). Strong gradients in DOM 
concentration exist between estuarine, shelf and open ocean waters (Hopkinson et al., 2002; Vlahos 
et al., 2002). The DOM concentrations on the shelf are elevated compared to the open ocean and 
contain a larger labile fraction and younger DOC (Bauer, 2002) than DOM in deep slope waters, 
where carbon is more refractory and enriched relative to nitrogen and phosphorus (Hopkinson et al., 
2002). Because half-lives of the labile DOM pool are on the order of shelf-water residence time, a 
substantial pool of DOM that is depleted in nitrogen and phosphorus relative to carbon remains and a 
net export of DOC to the open ocean can occur by advective and eddy diffusive processes 
(Hopkinson et al., 2002). A DOC budget study based on field measurements estimated a total export 
from the MAB shelf to the open ocean of 18.7-19.6 Tg.C.yr' 1 (Vlahos et al., 2002). To our 
knowledge, there has been no attempt to model the DOC dynamics and fluxes in the MAB. 
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DOM is believed to play an important role in carbon and nitrogen cycling from regional to global 
scales. DOM production is known to confound eutrophication issues (e.g. mucilage events) and may 
contribute to climate regulation by impacting carbon cycling and sequestration. Models have 
included DOM to study eutrophicated (Lancelot et al., 2005), eutrophic (upwelling, Ianson and 
Allen, 2002), mesotrophic (Anderson and Williams, 1998; Fasham et ah, 1999) and oligotrophic 
(Anderson and Pondaven, 2003; Raick et ah, 2005) regional systems and the global ocean (Popova 
and Anderson, 2002). Fasham et ah (1999) demonstrated that accounting for DOC is essential for 
euphotic ecosystem models and development of a carbon budget. Because DOC and DON dynamics 
are partially decoupled (e.g. carbohydrate production, mineralization rates), these models explained 
important deviations from the Redfield ratio in terms of productivity and export fluxes, and 
sustained nutrient-based primary production through DON mineralization and atmospheric inputs 
(Seitzinger and Sanders, 1999). 

The goal here is to describe the main pathways of DOM from production to mineralization in the 
MAB and to estimate the horizontal export of semi-labile DOC to the open ocean. This paper is 
organized as follows. After a brief description of the coupled physical-biogeochemical model, the 
DOM modeling is detailed (Section 2), the model validation and results will be presented (Section 3) 
and discussed (Section 4) with a focus on DOM dynamics and carbon export to the open ocean. This 
work is part of the Interdisciplinary Project USECoS, Eastern U.S. Continental Shelf Carbon 
Budget, for which a complete overview is presented in Hofmann et al. (2008). 

2. Model description 
2.1 The physical model 

The three-dimensional ocean circulation model (ROMS, Regional Ocean Modeling System 
version 3) extends across the Northeast North American (NENA) shelf including the Scotian shelf, 
the Gulf of Maine, the MAB, the South Atlantic Bight (SAB), the Eastern Gulf of Mexico and the 
adjacent deep sea (Figure 1). The focus of this paper is on the MAB and Georges Bank regions of 
the continental margin (Figure 1). The NENA model is nested within a North Atlantic (NA) model 
in order to capture large circulation features and variability such as the Labrador Current, the Gulf 
Stream and associated subtropical gyre circulation. 

ROMS (Haidvogel et al., 2000, Haidvogel, 2007) is a model widely used for shelf circulation and 
coupled physical-biological applications (e.g., Dinniman et al., 2003; Marchesiello et al., 2003; Peliz 
et al., 2003; Fennel et al. 2006; Wilkin, 2006). The ROMS computational kernel produces accurate 
evolution of tracer fields — a particularly attractive feature for biogeochemical modeling because it 
facilitates accurate interaction among tracers and accounting of total nutrient and carbon budgets. 
The application on the NENA domain uses a 10-km horizontal resolution and 30 terrain- following 
vertical levels stretched to give high resolution in surface and bottom boundary layers. This 
resolution is sufficient to capture the dominant dynamics governing shelf-wide circulation, yet is 
coarse enough to allow multiple simulations that explore sensitivities within the ecosystem model. 
Open boundary temperature, salinity and sub-tidal frequency velocity are taken from 3 -day averages 
of the Hybrid Coordinate Ocean Model data assimilation product developed as part of the Global 
Ocean Data Assimilation Experiment North Atlantic Basin ‘best-estimate’ analysis for 2003-present. 
Tides have been introduced at the boundary using harmonic data from the Oregon State University 
Topex/Jason altimeter data inversion and a surface gravity wave radiation scheme (Flather, 1976). 
Air-sea heat and momentum fluxes are computed using bulk formulae (Fairall et al., 2003) applied to 
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model sea surface conditions, and 
air temperature, pressure, humidity 
and winds from daily average 
National Center for Environmental 
Prediction re-analysis fields. 
Vertical turbulent mixing closure 
uses the parameterization of 
Mellor and Yamada (1982) and 
Warner et al. (2005). Coastal 
freshwater inputs are applied using 
USGS river flow data. 

This model exhibits recognized 
features of local and remotely 
forced circulation on the shelf and 
slope. These include wind-driven 
upwelling in the MAB and SAB, 
buoyancy-driven river plumes, 
tidal mixing and tidal residual 
mean currents on Georges Bank, 
southwestward mean flow in the 
MAB, low salinity on the MAB 
inner shelf, and retention of 
passive particles in the shelf-slope 
front. In addition the model 
captures Gulf Stream intrusions in 
the SAB, and interactions of Gulf 
Stream warm rings with the New 
England slope (Hofmann et al. 
2008). 

The simulations described below 
use a higher background value for 
diffusivity (10° m 2 s' 1 ) than in 
Fennel et al. (2006, 10' 6 m 2 s' 1 ) to 
compensate for the lack of 
secondary mixing processes such 
as internal waves and sub-grid 
turbulent diffusion. Similar 
background levels of diffusivity 
were shown to be necessary to 
reproduce the vertical temperature 
field and have been used in other 
model applications as well (10' 
~ m s’ , Kantha and Clayson, 

1994, and 2.1 0' 5 m 2 s' 1 , Ledwell et al., 1993). 


NENA domain 



Figure 1 . Geographical location and bathymetry (in meters) of the NENA 
(upper graph) and details of the MAB and Georges Bank regions (bottom 
The three numbered stations cited in the text are also shown. The locations 
the text are: ‘CC’ Cape Cod, ‘Co’ Comiecticut River, ‘Ho’ Housatonic Ri\| 
Long Island, ‘Hu’ Hudson River, ‘De’ Delaware Bay, ‘Ch’ Chesapeake B 
Roanoke River. 
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2.2 The biogeochemical model with DOM 

A general overview of the biogeochemical model schematic is presented in Figure 2 and the full 
details of the governing equations including the semi-labile DON and DOC are presented in the 
Appendix. The biogeochemical model includes the dynamics of nitrogen and carbon cycling. The 
semi-labile DOC and DON were added as state variables to the nitrogen-based model of Fennel et al. 
(2006) and the inorganic carbon component described in Fennel et al. (submitted). The nitrogen 
model includes eight state variables: nitrate (. NO 3 ), ammonium ( NH4 ), phytoplankton ( Phy ), semi- 
labile DON (named hereafter DON), zooplankton (Zoo), small and large detritus (SDetN and LDetN) 
and phytoplankton chlorophyll (Chi). The nitrogen-cycling formulations used here are the same as 
those in Fennel et al. (2006), except for several processes and parameters such as the resuspension 
and burial of POM (see Appendix and Table A2) which were required to meet the requirements of 
the DOM equations. The carbon model includes dissolved inorganic carbon (DIC), semi-labile DOC 
(named hereafter DOC) and small and large detritus (SDetC and LDetC). Phytoplankton and 
zooplankton are expressed implicitly in carbon units using the nitrogen unit equation and their 
specific C to N ratio (CNp and CNz respectively), thus no explicit equations are required (see 
Appendix). The DIC dynamics and air-sea exchange of carbon dioxide are described in Fennel et al. 
(submitted). 
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Figure 2. Flow diagram of the ecosystem model including the carbon (solid) and nitrogen (dash) cycles. 
The DOM sink and source fluxes are highlighted in red. 
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Semi-labile DOC and DON lability 

The definition of the DOC and DON pools may vary significantly between authors, therefore a 
definition is provided here in relation with the scales involved and goal pursued. The DOM pool is 
generally divided into labile, semi-labile and refractory pools although a continuum of biological 
lability exists between these categories (Amon and Benner, 1996). Highly variable decomposition 
(or turnover) rates of DOC were measured for surface and bottom waters of the MAB (Hopkinson et 
al., 2002). The refractory pool has a very long turnover time (several thousand years on average, 
Druffel et al., 1992; Santschi et al., 1995) and its concentration is relatively constant in the surface 
ocean at the yearly time scale. Although the refractory DOM represents 71% of the total DOC pool 
and 61% of the total DON pool in the shelf and slope water of the MAB (Hopkinson et al., 2002), 
the model does not take into account its variability since this study concentrates on the seasonal 
production of DOM. The labile material is defined here as having a turnover time scale of a few days 
to hours. Since it is mineralized in a few days within the 10 km-grid box of the model, the labile 
DOM is directed to the dissolved inorganic compartments (DIC and NH 4 ). The semi-labile fraction 
simulated by the model has a turnover time of one week to several months (due to a temperature 
dependency), which is on the order of the shelf residence time in the MAB (~ 100 days). As such 
this defined semi-labile DOM can therefore be efficiently exported to the open ocean by horizontal 
transport. 

DOM production by phytoplankton 

An overview of the literature highlights two phases of DOC production by phytoplankton. 
Sondergaard et al. (2000) suggest in their study that exponentially growing communities produce the 
most labile DOC, whereas declining and nitrogen-deficient communities produce the least labile 
DOC. During the growth phase, DOM production is linked to biomass and dominated by the 
exudation of labile-low molecular weight (LMW) organic compounds (Biddanda and Benner, 1997; 
Jensen, 1983; Lancelot, 1984) with a C to N ratio of ~7 (range 3-1 1 depending on species [—6.6 for 
the diatom sp. skeletonema ], Biddanda and Benner, 1997). In fact, the exudation was shown to be a 
passive diffusion across the outer cell membrane that occurs as long as new products of 
photosynthesis are available (Maranon et al., 2004). During the stationary and decaying phase of the 
bloom (i.e. under nutrient stress), large quantities of semi-labile, high molecular weight (HMW) 
DOM with high C to N ratios (10-25 compiling results from Benner et al. (1992)) would be released 
as a result of the exudation of polymeric carbohydrates (Lancelot and Billen, 1985) or due to cell 
lysis and to ‘sloppy’ feeding by zooplankton. 

Some evidence suggests that the release of carbohydrates by phytoplankton could mainly explain the 
accumulation of semi-labile DOC after the spring bloom and its progressive remineralization during 
summer and autumn. Biddanda and Benner (1997) showed that the relative abundance of 
carbohydrates in phytoplankton DOC increased from 23% during the exponential phase to 80% 
during the decay phase. Continued maintenance of photosynthetic machinery after nutrient 
exhaustion was found to be accompanied by excretion of DOM and especially carbohydrates with 
high C to N ratio (Hellebust, 1965; Norrman et al., 1995). Diatoms can continue to excrete 
polysaccharides for a considerable time after the halt of cellular protein synthesis (Jensen, 1983). In 
many offshore systems, a DOC decrease is found to continue after nutrient exhaustion (Sambrotto et 
al., 1993). The semi-labile DOC release would thus occur mainly in low nutrient conditions and is 
likely to be associated with phytoplankton primary production. Furthermore, Thingstad et al. (1997) 
proposed that the spring to summer accumulation of DOC is related to microzooplankton grazing on 
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bacteria coupled to low bacterial growth rates, which would reduce DOC remineralization and allow 
DOC accumulation. 

Two models of extracellular DOM release have been proposed: the overflow model (Fogg, 1966; 
1983; Nagata, 2000; Williams, 1990) and the passive diffusion model (Bjomsen, 1988; Fogg, 1966) . 
Even if these models were opposed in conflicting reports, it is likely that they are not mutually 
exclusive and that both models are correct given the right environmental conditions and plankton 
community structure (Carlson, 2002). Carlson suggests that the extracellular release of labile-LMW- 
DOM model is likely to be a passive diffusion process linked to biomass. During the stationary and 
decay phases, the overflow model is likely to represent an active release of semi-labile-HMW-DOM 
linked to primary production and enhanced in a nutrient-depleted environment. Because primary 
production is traditionally expressed in models as a function of biomass, both terms of exudation 
(labile) and excretion (semi-labile) of DOM are dependent on primary production in the present 
model (and in most other modeling studies, e.g. Anderson and Williams, 1998). This option allows 
in addition an easier interpretation of the model results. 

Exudation of labile DON and nutrient-based labile DOC 

Admiraal et al. (1986) estimated that approximately 3% of assimilated nitrate is excreted as amino 
acids. The rate chosen for the labile DON exudation (and instantaneous mineralization in the model) 
is set to con = 3 % of phytoplankton nitrogen production. The labile DOC leakage is also expressed 
as a function of primary production with the same rate ( co c = 3%) to ensure a constant Redfield ratio 
for labile DOM and phytoplankton. 

Semi-labile DON exudation by phytoplankton 

Bronk et al. (1994) found an average DON release of 25 to 41% of the inorganic nitrogen uptake in 
offshore oceanic (25%), coastal (27%) and estuarine (41%) environments with turnover times of 10 
±1,18 ±14, and 4 days respectively. Sondergaard et al. (2000) found 25% in their mesocosm study. 
Since the labile fraction of DON production is estimated to be a few percent of DIN uptake (see 
above), the semi-labile DON total release in the coastal ocean is most probably in the range of 22 to 
24% of nitrogen-based primary production with decreasing value towards offshore. Varela et al. 
(2003) provided some evidence that DON production is dominated by grazing processes rather than 
by direct phytoplankton excretion. Large DON losses (>50% of nitrogen uptake) were attributed to 
intense grazing and sloppy feeding for several marine ecosystems (Bronk and Ward, 2000). The 
maximum of DON release was found to occur when small, presumably hetero trophic, flagellates 
dominated the biomass and not the primary production (Varela et al., 2003), i.e. sloppy feeding by 
flagellates could significantly increase the DON release. In summary, sloppy feeding might 
dominate the DON release during a short period of intense grazing, but phytoplankton exudation and 
detritus solubilization dominates otherwise. It is estimated that the DON released by exudation 
follows a Redfield ratio of the nutrient based DOC exudation that is set to 4% of primary production 
(basal value of DOC exudation by healthy phytoplankton, see next subsection). The rate of semi- 
labile DON exudation by phytoplankton (ef) is set to 4% of nitrogen-based primary production. 
Following the above assumption that semi-labile DON release should be about 22% of primary 
production, the sloppy feeding should account for about 2 to 14% (low and high grazing) and PON 
solubilization for about 16 to 4% depending on grazing. 
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Nutrient and carbon excess-based semi-labile DOC exudation by phytoplankton 
In their mesocosm experiment, Norrman et al. (1995) observed that 23% of total new production 
accumulated as DOC, which was found to increase due to a combination of excretion and cell lysis. 
A large range of values of DOC production as a fraction of primary production can be found in the 
literature (5-30%, Biddanda and Benner 1997, Norman et al. 1995), however excretion from natural 
healthy phytoplankton was found to be lower (4-16%) than at the end of a diatom bloom (17-38%, 
Hellebust 1965). 

Two terms describe the semi-labile DOC exudation by phytoplankton in the model: a nutrient-based 
and carbon excess-based release. The nutrient-based release reflects the healthy phytoplankton 
exudation of semi-labile DOC and follows the semi-labile DON exudation with the Redfield ratio. 
The carbon excess-based release represents the carbohydrate over-production by nutrient- stressed 
cells. The carbon excess uptake is seen as an ‘overflow’ of photosynthesis under nutrient limitation. 
It is formulated as the difference between the nutrient-saturated (light-limited) and nutrient-limited 
(light-limited) primary production and is directed to the semi-labile DOC (Andersen and Williams 
1998, Ianson and Allen 2002, see Figure 3 with the details of the terms in Table 1). The carbon 
excess uptake (t/exc.c) is thus expressed: 

Uexc.c = 7 CNp ( PPl — PPl Ln ) 

where PP L is the nutrient-based primary production limited by light, L N is the nutrient limitation, 
CNp is the C to N ratio for phytoplankton and y the parameter of carbon excess-based DOC excretion 
by phytoplankton. A fraction (o c ) of the carbon excess uptake is directed to the semi-labile DOC 
pool and represents the exudation of carbon excess-based DOC release. This fraction is set to 
cr c =0.45 (Biddanda and Benner, 1997 found -35%). The labile DOC originating from carbon excess 
uptake represents a slightly higher fraction (l-o c ) and is directed back to DIC. 



The total excretion is commonly expressed as a fraction of the total carbon fixed by phytoplankton, 
the percentage extracellular release (PER). In the present setting, the PER follows: 
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The PER for diatoms was estimated to be between 10% and about 55% (Baines and Pace, 1991; 
Obemosterer and Hemdl, 1995) with an increase between the exponential and the stationary phase of 
the bloom. Higher PER values (70-80%) were observed in eutrophic water for Phaeocystis pouchetii 
(Lancelot, 1983). The analytical formulation of the PER in our model shows that the mean PER is 
lower than 65% for Ln below 0.5 (nutrient limiting condition) when y is set to 0.20. We thus chose y 
=0.20 for the simulations. For comparison, Anderson and Williams (1998) adjusted y to 0.26 to 
achieve the desired spring DOC concentration and obtained PER values between 1 0% and 60% for 
Station El in the English Channel. 

DOC release by ‘Sloppy’ feeding 

Measurements from the literature suggest a high DOC release when the prey is large relative to the 
copepod and low DOC release when the prey is small relative to the copepod (Moller, 2005). During 
a diatom bloom, sloppy feeding was, by far, the most important contributor to the DOC production 
by Calanus spp., and 49% of the carbon removed from suspension by the copepods was returned to 
the water column as DOC (Moller et al., 2003). Moller (2005) found a significant relationship 
between the DOC production through sloppy feeding by zooplankton and the copepod-to-prey size 
ratio. Q defines the fraction of prey carbon removed from suspension and lost as DOC for 
copepod:prey size ratio below 55: 
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where ESD is the equivalent spherical diameter. 


According to Moller et al. (2003), when copepods graze large diatom cells in spring the 
copepod:prey size ratio can reach a minimum of 10 and Q values may reach 71%. For an increasing 
size ratio, i.e. when copepods graze on smaller prey during summer, Q decreases linearly down to 
£9=0.1% for a size ratio of 55 (or more). In agreement, Moller (2005) illustrates that when the prey is 
large relative to the copepod, i.e. during a bloom of large cells, copepods lose significant amounts of 
dissolved material. In contrast, the link between copepod feeding and energy flow to higher trophic 
levels is tighter when the prey is small, i.e. during oligotrophic periods when small cells dominate 
the phytoplankton prey. The fraction of DOC released by sloppy feeding is likely to reach its 
maximum during the spring bloom (large diatoms) and minimum in summer when smaller cells are 
grazed. Since a high grazing level is a good proxy of high biomass of large cells (diatom spring 
bloom), a linear relationship is used to enhance the fraction of DOC release by sloppy feeding at 
high grazing levels: 
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where g is the grazing, g max is the maximum grazing rate, g = g max (Fennel et al., 2006) 

k p + Phy ‘ 

and k p is the half- saturation constant of phytoplankton ingestion. The maximum fraction of DOC 
release of 71% is thus encountered when the grazing intensity is maximum, i.e. at the highest levels 
of phytoplankton biomass. 

The fraction of semi-labile DOC (to total DOC) in the phytoplankton cell is estimated using the 
work of Biddanda and Benner (1997). They estimated that dissolved carbohydrates represent the 
major part of cell DOC during the stationary phase of the bloom. They measured a fairly stable 
fraction of dissolved polysaccharide carbohydrate (between 78 and 94% with a mean of 85%) 
compared to dissolved monosaccharide carbohydrate between the exponential growth and the 
stationary phase. However, dissolved polysaccharides (particularly fresh material) can be rather 
labile (Mannino, 2000; Mannino and Harvey, 2000a), therefore the fraction of semi-labile DOC 
compared to labile DOC is set to 5c=55%. Although the cell-content of DOC is constant in the 
model, Biddanda and Benner (1997) measured an increase of dissolved carbohydrates from 23 to 
80% of cell DOC for four phytoplankton groups ( Synechococcus , Phaeocystis , Emiliana and 
Skeletonema ) from the growth to the decay phase. The sloppy feeding related terms for carbon are 
therefore the following: 

Assimilation of organic carbon by zooplankton 
Semi-labile DOC release by sloppy feeding 
Labile DOC (towards DIC) release by sloppy feeding 
Fecal pellets production (to large detritus) 
where (3 is the zooplankton assimilation rate. 

DON release by ‘sloppy’ feeding 
Hasegawa et al. (2001) found that 9 to 75% of ingested nitrogen is assimilated in zooplankton 
biomass depending on food concentration. Therefore, from 25 to 91% of grazed nitrogen is released 
as PON, DON or ammonium for low and high food concentrations respectively through the 
processes of sloppy feeding, excretion and egestion of fecal pellets. Zooplankton excretion rates in 
the original model (Fennel et al., 2006) included the assimilation related excretion (/#= 0.1 d" 1 if 
grazing is maximum) and the basal metabolism related excretion (1bm = 0.1 d’ 1 ). The rate of fecal 
pellet production was set to 25% and the assimilation efficiency to 75%. Therefore, in order to 
compare with the values of Hasegawa et al. (2001), the fraction of nitrogen released per day in the 
model by excretion (assuming zooplankton ingests 60% of its weight of prey per day 1 ) and egestion 
of fecal pellets was 0.1*0.75+0.1*0.6+0.25=0.385. This value can decrease to 28.5% if the 
zooplankton stops grazing. Therefore the model absolute assimilation efficiency ranged from 61.5% 
for high food condition to 71.5% under low food condition, which is a low range for DON release 
compared to Hasegawa et al. (2001, 9-75% assimilation). Adding the process of sloppy feeding has 
the effect of decreasing the absolute assimilation efficiency for high food condition, with a 
maximum contribution of about 50% of the nitrogen grazed. The assimilation efficiency in the 
model was constant and linked to fecal pellet production (/?=0.75). Compared to the original model, 
p is grazing-dependant taking into account the loss of DON by sloppy feeding. IN the current model, 
the ‘absolute’ assimilation efficiency (which includes the excretion rates, p - 1 e P g/gmax - Ibm ) varies 
between 10 and 75% as a function of the grazing intensity (Hasegawa et al., 2001), which leads to: 


= CNp p g 

= CNp (1-p) Qdoc be g 
= CNp (1-p) Qdoc ( 1-dc) g 
= CNp (1-p) (1 -Qdoc) g 


1 The maximum grazing rate is 0.6 d+ 
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where fS is the grazing-based assimilation efficiency (excluding the excretion rates). The assimilation 
efficiency defined here does not include the excretion rates and varies between /?= 0.22 for g=g,„ax 
and /?=0.85 for g=0. The fraction of semi-labile DON to total DON is assumed to be low (the 
opposite for DOC) and is set to (5y=30%. Qdon being the fraction of total DON to (DON+PON) 
within the phytoplankton cell, a fraction (1 -Qdon) of the remaining non-assimilated material (l-fi) is 
allocated to the fecal pellets and the complementary fraction ( Qdon ) lost by sloppy feeding to the 
DON. Seventy percent (1 -S N ) of this last fraction is labile and 30% ( =S N ) is the semi-labile DON. 

In summary, the grazing term for nitrogen is divided as following: 


Nitrogen assimilation for zooplankton 
Seini-labile DON release by sloppy feeding 
Labile DON release by sloppy feeding 
Fecal pellets production (to large detritus) 


= Pg 

= ( 1 -/0 Qdon Sn g 
= ( 1 ~P) Qdon ( 1 -Sn) g 
= (!-/?)(! -Q D ON)g 


The C to N ratio of semi-labile DOM produced by sloppy feeding in the model, deduced analytically 
using the terms defined above, has the constant value of 12.1 (CNp S c / S N ) using the current 
parameter set (see Appendix). 


Solubilization of carbon and nitrogen detritus 

Smith et al. (1992) showed that bacteria attached to aggregates can express high levels of hydrolytic 
ectoenzymes which results in the release of DOM. They suggest a preferential solubilization of 
nitrogen detritus over carbon detritus in agreement with the observed increase in C to N and C to P 
ratios of sedimenting material with depth. In the experiments, the aggregates released mainly 
dissolved combined amino acids (DCAA). Solubilization was uncoupled from bacterial uptake, with 
50-98% of the DCAA released escaping rapid utilization from attached bacteria. From the data 
estimated by Smith et al. (1992, Table 2), DON release rates from aggregates are calculated 
assuming a C to N ratio for POM of ~5 as observed by Harvey et al., 1995). The solubilization rates 
for DON range from 0.066 to 0.200 d' 1 with a mean of 0.12 d" 1 .. The solubilization rate of the large 
and small carbon detritus is set to ssDetc = SLDetc = 0.08 d’ 1 . The higher rate chosen for the nitrogen 
detritus solubilization is SLDetN = ssDetN = 0.11 d’ 1 consistent with the carbon detritus specific 
respiration rate on diatom aggregates of 0.083±0.034 d' 1 (Ploug and Grossart, 2000). The fractions of 
semi-labile DON and DOC cell-contents (Sn and Sc) described above are used to quantify the release 
of semi-labile DOM by solubilization. The labile DOM generated by solubilization (directed to the 
ammonium and DIC pools) uses the complementary fractions (1 -Sn) and (l-<Sc). 


Mineralization of semi-labile DOC and DON 

Since the photo-oxidation process is likely to impact refractory DOC, it is not included in the model. 
Hopkinson et al. (2002) estimated rates for DOC and DON mineralization by bacteria at 19-20°C in 
the MAB. The mean DON degradation rate measured (see Table 4 in Hopkinson et al. 2002) with 
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half-life values comprised between 4 days (0.173 d' 1 ) and 365 days (0.0019 d' 1 ) is 0.073±0.043 d' 1 
with a corresponding value for DOC of 0.055±0.057 d' 1 . In a recent six-month degradation 
experiment at 19-20°C with DOM collected in the MAB, rates of DOC mineralization are between 
0.015 d' 1 and 0.043 d' 1 (Russ and Mannino, in prep). The values proposed here are 0.029 d' 1 (ti /2 = 
24 d) for the semi-labile DOC and 0.060 d' 1 (ti /2 = 12 d) for the semi-labile DON at 19-20°C. Since 
temperature dependence is likely to occur with a Qio = 2 (Pomeroy et al., 1991), the degradation rate 
has the following formulation: 

Mineralization rate [d' 1 ] = ao e 0 07T 

where ao= accT 0.00767 d" 1 for carbon, acr um? = 0.0153 d" 1 for nitrogen semi-labile DOM and T is the 
temperature in °C. Consequently for the semi-labile DOC, the degradation rate ranges between 
0.0088 d' 1 at 2°C to 0.0545 d' 1 at 28°C and between 0.01 76 d' 1 at 2°C to 0.1 086 d' 1 at 28°C for the 
semi-labile DON. Although bacteria are involved in DOM degradation, bacteria are not explicitly 
included in the model to avoid the multiplicity of parameters that cannot be calibrated using our 
dataset. However, bacterial processes such as organic matter mineralization and solubilization of 
POM are included in the model with a Qio temperature formulation. 

2.3 Other modifications to the model 

An additional light attenuation coefficient to account for colored dissolved organic matter (CDOM) 
was added to the water and chlorophyll attenuation in response to an overestimated euphotic depth 
simulated by the model in the MAB compared to observations. The light attenuation formulation for 
CDOM absorption correlated to salinity was implemented in agreement with the combined CDOM 
plus detritus absorption measurements at 443 nm (acDM(443), Magnuson et al., 2004) in estuaries 
(Chesapeake Bay: 0.368 ± 0.076 m" 1 in the mid-Bay for a salinity of 12.9 ±3.2 psu, 
0.284 ± 0.090 in' 1 in the lower Bay for a salinity of 20.8 ± 4.9 psu), on the shelf (0.168 ± 0.057 m' 1 
in the inshore MAB and 0.070 ± 0.035 m" 1 in the offshore MAB) and in oligotrophic waters 
(0.010 d" 1 at the BATS station). The formulation used is (with salinity in psu): 

Kcdom = 0. 5329-0. 02669 *salinity+0. 0003395 *( salinity ) 2 

It is slightly lower (~0.1 m' 1 for salinity lower than 30 psu) than the observed acDM(443) to account 
for a lower absorption of the PAR in the entire visible spectrum than in the blue (443 nm). 

In the present model version, the zooplankton products (dead zooplankton and fecal pellets) and the 
aggregates of dead phytoplankton cells with small detritus are directed to the large detritus 
compartments instead of the slow-sinking small detritus compartments in Fennel et al. (2006). As a 
result in the water column, the small detritus compartments ( SDetN and SDetO) are only fed by the 
dead, non-aggregated, phytoplankton cells and the large detritus compartments (LDetN and LDetC ) 
include the aggregates of dead phytoplankton cells with the small detritus and the zooplankton 
products. 

In order to explore the high POM burial of the U.S. northeastern continental shelf (Thomas et al., 
2002) and to compare it with the horizontal export of DOM to the open ocean, a parameterization of 
POM resuspension and burial was added to the model’s bottom boundary formulation. The 
resuspension rate of the POM flux reaching the seabed is a function of the bottom friction velocity. 
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The resuspended fraction of POM is thus largely dependent on the local near bottom current velocity 
that is driven by the general circulation and tides on the continental shelf and by wind events in 
shallow waters. The remaining fraction of PON and POC accumulates and is buried in the sediment 
assuming that the burial efficiency of the particulate organic carbon is proportional to the vertical 
flux of POC reaching the seabed (Henrichs and Reeburgh, 1987). The POM resuspension and burial 
processes are fully described in the Appendix. 

2.4 Initial and boundary conditions 

Initial and boundary conditions for nitrate were derived using polynomial approximations that 
predict nitrate concentration from temperature using the NODC World Ocean Database 2001 
(Fennel et al., 2006). The semi-labile DOC and DON boundary conditions are constant and set to 

1.0 mmol m" and 0.15 mmol m" respectively. This approximation does not affect the MAB area 
because it is far enough from the boundary limits of NENA (Figure 1) to lose its memory (the time 
required to transport water masses from the boundary limit to the MAB is lower than the semi-labile 
DOM halftime). The same reasoning was applied to all other biological state variables with 
boundary conditions set to small background values. Monthly climatology for river flow, nitrate, and 
ammonium were derived from the U.S. Geological Survey monitoring database. The MAB receives 
a large supply of freshwater from the Hudson, Delaware and Chesapeake estuaries that carries high 
loads of DOC. Much of this DOC is believed to be of terrestrial origin, consisting of mostly 
refractory organic matter that can be transported across the shelf and into the open ocean (Aluwihare 
et al., 2002; Bauer, 2002). Although the semi-labile fraction is believed to be small, the main input 
of freshly produced DOM to the MAB is from production within the estuaries. Since the boundary 
conditions for the Chesapeake Bay and Delaware Bay are located in the mid and upper bay 
respectively, a polynomial fit is applied to the DOC data available in these areas from the 
Chesapeake Bay Program database (http://www.chesapeakebay.net/wquality.htm) and other sources 
(Jon Sharp, pers. comm.). A dampened seasonal pattern is emphasized with a minimum DOC value 
in March and a maximum value in early October with a mean of 297±14 mmol m" . A small fraction 
(10%) represents the semi-labile DOC and is used as a boundary condition for rivers. A Redfield C 
to N ratio, which approximately characterizes freshly produced DOM (i.e. the semi-labile fraction), 
is used to derive the semi-labile DON boundary condition for rivers. The river boundary condition 
for phytoplankton biomass and chlorophyll is set to background values of 1.8 mmol N m' and 

6.0 mg Chi m’ respectively. 

3. Model results 

A spin up of seven months is conducted on the NENA domain to initialize the biogeochemical 
model (June 2003) for the year 2004. The results of twin simulations -with and without the DOM - 
are discussed in the following subsections. 

3.1 Vertical and seasonal distribution of DOM 

A station chosen for the representation of the vertical and seasonal DOM processes from the 
southern MAB shelf is detailed in this section (station 3, see location in Figure 1) and compared to 
results in the central MAB and Georges Bank areas. As revealed by the simulation, station 3 is 
located in a region of marine production and export of DOC. The water depth at that station is only 
46 m, but this area is located near the shelf break at the vicinity of the Gulf Stream. The annual 
primary production at this site, 25 1 g C in'" yr' , is characterized by a nitrate-sustained production at 
the subsurface in summer (Figure 5a and 5b). A nitrate concentration of 5 mmol m' 3 is encountered 
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at the depth of the 10% isolume during that period following an intrusion of the Gulf Stream as 
suggested by Schollaert et al. (2003) and observed in July 1996 (Redalje et al., 2002). The most 

•31 T 

important period of nutrient uptake, with values up to 0.8 mmol NO3 m‘ d' and 0.5 mmol NH 4 m' 
d" 1 (Figure 5a), occurs during summer between the mixed layer depth and the depth of the 10% 
isolume. The annual new production is 46% of the nutrient-based production (nitrate and 
ammonium), compared to 36% at station 2 where the Gulf Stream has less influence and 41% at the 
well -mixed tidal-driven station 1. The carbon excess uptake at station 3 of 30 g C m‘ yr" represents 
12% of total primary production (14% at station 2 and 6% at the nutrient-rich station 1). This ‘extra’ 
carbon uptake occurs at a shallower depth than the nitrogen-based uptake, within the nutrient- 
depleted mixed layer and immediately below where the gradients of biomass and light intersect. The 
spring bloom arises in March and early April at station 3 (Figure 1 la), in late March and April at 
station 2 and in May at station 3. Contrary to a traditional maximum of phytoplankton biomass in the 
mixed layer during the spring bloom (station 2), the maximum biomass at station 3 is located at the 
subsurface during summer in relation with the Gulf Stream intrusion. 

The bulk of semi-labile DON is formed at the three stations during the bloom and post-bloom 
periods with a maximal content near the surface of 3 mmol m" . This value is in agreement with the 
seasonal increase observed by Hopkinson et al. (1997), from 2 to 5 mmol m‘ ). The semi-labile DOC 
pool at station 3 increases rapidly during the decay phase of the spring bloom with values up to 
55 mmol m"' (Figure lib). This freshly produced DOC pool is sustained until the end of May in 
relation with a high subsurface phytoplankton biomass. The progressive decrease of the subsurface 
biomass in June and July leads to a decrease of the surface semi-labile DOC content from 40 to 
15 mmol m' 3 . A second increase to 23 mmol m° is observed in August linked to a shallower 
stratification. A slight subsurface maximum of semi-labile DOC and DON is observed at station 2 in 
July and August (results not shown). Homogeneous profiles of DOM are found at station 1 , which 
alternate with surface maxima in relation with periodic stratification events. Figure 10 shows the 
seasonal variability of the semi-labile DOC concentration at the surface in the MAB and Georges 
Bank. The highest values in the MAB occur in spring (-65 mmol m‘‘) and in summer in Georges 
Bank with highest value in the shallowest area (from 30 to 50 mmol m" ). The semi-labile DOC 
concentrations remain elevated in rivers, estuaries and plumes of the MAB with, however, lower 
values in autumn. The slope off the MAB and Georges Bank between 100 and 1000 m shows a local 
maximum concentration in spring correlated with high primary production levels. 

The model exhibits a stable and lowest C to N ratio of DOM during the growing phase of the spring 
bloom and an increasing ratio during the stationary and decaying phase of the bloom. The increase 
of C to N ratio (atoms) of DOM during the post-bloom in spring is from 10.5 to 14. 5 at station 1, 12.5 
to 21.5 at station 2 and from 12.5 to 18.5 at station 3 in agreement with previous measurements 
(between 10 and 25, Benner et al., 1992). 

3.2 Horizontal DOM distribution and model evaluation 

The general distribution of total DOC concentration in surface waters described by Vlahos et al. 
(2002) shows an increase from northeast to southwest, and from offshore to inshore. These gradients 
arise from the production and accumulation of total DOC concentration on the shelf as a result of 
primary production and river inputs to southwestward flowing water mass (Vlahos et al., 2002; 
Mannino et al. 2008). The validation of the modeled semi-labile DOC is complicated by the different 
quantity measured in the field which includes also the labile and refractory fractions. Furthermore 
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due to the particularly coarse model resolution in the estuary areas, we have less confidence in the 
simulated estuarine production, which is mostly exported in the inner- and part of the mid-shelf. We 
focus here on the outer shelf and slope where the model is well adapted and the exchange with the 
open ocean is high. Vlahos et al. (2002) derived areas of DOC production in the MAB subtracting 
from the total measured DOC the “expected” DOC concentration based on conservative mixing 
(using salinity). The “excess” DOC derived outside of the river plumes, i.e. mainly in the outer-shelf 
and slope, is thus comparable to the simulated semi-labile DOC. Vlahos et al. (2002) found an 
excess DOC from 20 to 55 mmol C m° in the surface outer and slope waters of the southern MAB in 
March 1996 and April 1994. In August 1996, values range from 10-15 to 25-30 mmol C m" 3 in the 
same area. The DOC decomposition study of Hopkinson et al. (2002) for near-surface waters of the 
MAB demonstrated a degradable DOC pool of 18-40 mmol C m° in March 1996 and 23- 
50 mmol C in"' in August 1996 with two DOC components with half-lives of 1-14 days and 8-147 
days. Mannino et al. (2008) attributed a seasonal increase in DOC of 12 to 34 mmol C m’ 3 within the 
southern MAB between Spring and Summer of 2005 and 2006 to net ecosystem production 
(phytoplankton exudation, particle solubilization, grazing, etc.) of semi-labile DOC. The model 
ranges from 15 to 60 mmol C m' in March-April and from 15 to 35 mmol C m"' in August. At 
Georges Bank, the DOC dynamics show no major accumulation in the surface water during spring 
(Figure 10) in agreement with the observation (Chen, 1996). 

Figure 4 presents the comparison of the mean June 2004 surface chlorophyll concentration derived 
by the model and by the satellite sensor SeaWiFS in the MAB and Georges Bank regions. The 
simulated chlorophyll distribution agrees well with the satellite observations particularly at the 
enhanced concentration on Georges Bank due to the tidal mixing and permanent nutrient availability. 
The river plumes and inner shelves are also well reproduced, including in the Chesapeake Bay where 
Harding et al. (2005) showed that SeaWiFS overestimates in situ chlorophyll by approximately 
100%. A feature which is not observed by the satellite sensor but appears in the model and is well 
documented (Ryan et al., 1999) is the enhancement of surface chlorophyll at the shelf break of the 
MAB and southern Georges Bank from mid-April to late June. These higher chlorophyll 
concentrations corresponds to the transition period from well-mixed to stratified conditions and is 
sustained by the upwelled nitrate-rich waters of the geostrophic jet that flows along the shelf break 
and slope from Georges Bank to Cape Hatteras (Ryan et al., 1999; Figure 5a and 5e). 

A novel plot for quantitatively evaluating and displaying the skill of coupled biological-physical 
models, called the target diagram, has been recently introduced (Joliff et al., 2007). In these 
diagrams, bias and centered-pattem RMS are normalized by the standard deviation of the 
observations and plotted on the x- and y-axes, respectively. Because the sum of the squares of these 
two components of the RMS difference is equal to the square of the total RMS difference, the 
distance from the origin to each plot symbol represents total RMS error. Although centered-pattem 
RMS is inherently a positive quantity, in the target diagram the centered-pattem RMS is multiplied 
by the sign of the difference: standard deviation of observations - standard deviation of model. Thus 
symbols are plotted with positive x-coordinates if the model overestimates the variability of the data, 
and with negative x-coordinates if the model underestimates the variability of the data. The circle 
representing total RMS difference =1.0 (i.e. total RMS equals the standard deviation of the 
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Figure 4. Surface chlorophyll concentration (mg Chi in' 3 ) simulated by the model (left) and observed by 
the satellite sensor SeaWiFS (right) using the algorithm OC4V4 for June 2004 in the MAB and Georges 
Bank regions. Isobaths are in meters. 


observations) is typically superimposed on these diagrams for reference. By definition, model results 
falling within this circle reproduce observed quantities better than the mean of those observations. 
Target diagrams for both model runs (with and without DOM) are presented in Figure 6, and 
illustrate the skill of the models in reproducing the satellite monthly (2004) mean surface chlorophyll 
concentrations for the region encompassing the mid- and outer-shelf and slope of the MAB. This 
figure shows that the reference run generally (except February -April) underestimates the surface 
chlorophyll content, whereas the DOM run more often overestimates surface chlorophyll, especially 
in the spring (March and April). In the summer (June - August) the bias is nearly zero for the run 
with DOM., whereas the model results are negatively biased in the run without DOM. In these 
summer months both models underestimate the observed variability to the same degree. Both models 
also similarly overestimate the spatial variability of surface chlorophyll in February - May. 

The organic carbon production in the MAB and Georges Bank regions (Figures 7c and 5c) range 
between 100 and 300 g C m" yr in general agreement with previously published ranges (Falkowski 
et al. 1988 and Berger 1989: 120-300 g C m"“ yr ) with little alongshore variability on the central 
MAB (O’Reilly and Busch, 1984). However, this level of productivity is lower than the approach 
using mixed satellite and in situ profiles used by Mouw and Yoder (2005) with values of 320 g C m" 
yr on the shelf, 304 g C m" yr on the shelf break and 411 gC m" yr on the slope of the 
northern MAB. Similarly higher productivity was estimated by O’Reilly et al. (1987) on Georges 
Bank with 455 g C m’ yr on the shallowest sector, 310 g C m" yr between the 60 and 100 m 
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isobaths and 265 g C m"“ yr between the 100 and 200 m isobaths. Compared to these higher 
estimates, the model estimates are too low by approximately 1 00 g C m' yr . 

The seasonal variability of simulated primary production on Georges Bank (Figure 8a) in May to 

September are in agreement with measurements (slightly above 1 g C m'~ d , O’Reilly et al., 1987), 

2 —1 

but the December to March levels of productivity of 0.01 to 0. 10 g C m"“ d in the model are 
significantly lower than observations (0.2 to 0.6 gCm' d ). The lack of a low-light sensitive 
phytoplankton group in the model is believed to be the cause of the winter underestimation of 
productivity. The percentage of extracellular release (PER) shown on Figure 8 is in agreement with 
the measurements of O’Reilly et al. (1987) between 14% in the shallow water (8 to 23% at station 1) 
and 2 1 % over the slope of Georges Bank. Daily primary production levels simulated at station 2 
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Figure 5. Seasonal profiles of (a) new, (b) regenerated production (mmol N m‘ 3 d -1 ) and (c) total carbon 
production (mmol C m' 3 d _1 ) simulated by the model at station 3 (southern mid-shelf MAB). The mixed 
layer depth (solid black) and the depth of the 10% (dash white) and 1% (solid white) isolume are also 
presented. 


station 3 (depth = 46 m) 


(Figure 8b) are also in agreement with the measurement of the SEEP-I experiment in 1984 with 
0.60 g C m~ 2 d _1 in March and 1.33 g C m" 2 d -1 in April (Falkowski et al., 1988). The inner shelf 
north and south of the Delaware Bay show a much stronger underestimation of productivity in the 
model (100 g C m’ 2 yr” 1 ) than the observation (505 g C m" 2 yr” 1 , O’Reilly et al., 1987). In that 
particular area, the phytoplankton is nitrogen-limited and the lack of sediment erosion by waves and 
POM resuspension is believed to cause such discrepancy. The levels of productivity provided by the 
model in the southern MAB shelf between the Chesapeake Bay mouth and Cape Hatteras (0.5- 

1.0 g C m" d in March and July) are in the range of field data (0.5- 1.0 g C m' d in March and 
0. 5-2.0 g C m' 2 d” 1 in July, Verity et al., 2002). With the exception of the inner-shelf of the MAB, 
the productivity is thus well reproduced for the spring-summer-autumn period when the main DOC 
production events occur. POC deposition at the seabed was fairly well represented in the MAB shelf 
and slope. The POC deposition rates on the shelf off Delaware Bay and Chesapeake Bay are 2.7 and 

2.1 mol C m’ yr , respectively, (Biscaye et al., 1994) in agreement with the model which ranges 
from 1.0 to 2.5 and from 1.0 to 4.0 mol C m'“ yr , respectively (result not shown). The model, 
however, does not reproduce the POC deposition on the slope, with an underestimation factor of ~10 
compared to field measurements (4.6 to 13.1 mol C m'“ yr off Cape Hatteras and 1 to 2 mol C m‘ 
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yr on other slope areas of the MAB; Biscaye et al., 1994; Schaff et al., 1992; Thomas et al., 
2002). This suggests that the sinking velocity of the large detritus is too low or a third ‘very large’ 
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detritus pool with sinking velocities of 100 m d 1 (Walsh, 1994) is lacking in the model to reproduce 
the POC deposition and burial on the slope. 


However, the focus in the present study is on the 
shelf burial of POC in comparison with the 
horizontal DOC export. Contrary to the POC 
deposition, the POC burial on the shelf is not well 
documented. Although the burial is probably 
slightly underestimated, consistent with a moderate 
underestimation of the primary production, we 
believe the model estimate of POC burial (Figure 
9b) is realistic since the POC deposition rate agrees 
with field measurements and the simulated POC 
burial on the shelf is globally higher (from 1.5 to 
4. 1 mol C m" yr ) than the field estimates on the 
slope (1 to 2 mol C m’ 2 yr -1 ). 

3.3 Semi-labile DOC export to the open 
ocean and POC burial 

The annual mean of horizontal divergence of semi- 
labile DOC integrated over the water column 
(Figure 9a) shows specific areas of production and 
export (negative values) and areas of import 
(positive values) for both the shelf and the open 
ocean. Areas of high primary production are 
identified as regions of significant export of semi- 
labile DOC. The DOC release in Georges Bank and 
the shelf south of Cape Cod is mainly exported in 
the central MAB shelf. The DOC released on the 
southern outer-shelf and slope of the MAB and off 
Cape Hatteras is exported to the adjacent deeper 
ocean. Both areas of export and import show an 
annual flux on the order of 1 mol C m’ yr' . 

The maxima in POC burial (Figure 9b) occur on the 
inner-shelf south of the dominant simulated rivers 
and estuaries (from North to South: the 

Connecticut, Housatonic and Hudson Rivers, the 
Delaware and Chesapeake Bay, the Roanoke River, Figure 1). The productivity level and the water 
depth are the primary factors that determine the flux of organic matter near the seabed. The bottom 
friction then controls the POM deposition and the POM flux at the seabed, which regulate the burial 
through the burial efficiency. The river-influenced areas of the inner-shelf show a POC burial from 
1.5 up to 4 mol C m'^ yr" . Except for the region south of Cape Cod and Georges Bank where the 
tidal-induced bottom friction prevents deposition, the rest of the shelf shows decreasing values of 
POC burial from inshore to offshore, with a flux of about 0.5 to 1 mol C m’ yr" following the 40 m 
isobath. For water depth greater than 100 m in the MAB, the POC is entirely mineralized in the 
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Figure 6. Target diagram (representing the total RMS 

difference, see text for details) of the monthly mean 

surface chlorophyll concentration between the model 

estimate and derived from satellite (SeaWiFS sensor) 

for the reference run 

(upper) and the run including 

DOM (lower). The area included in this evaluation 

comprises the Mid- and Outer-shelf and the slope of 

the MAB. 



19 





300 


250 


200 


150 


100 


50 


yr 


-i 



1 300 


250 


200 


150 


100 


50 


76 -74 -72 -70 -68 -66 „ ^ -2 , -1 


g C m yr 


50 


40 


30 


20 


Figure 7. Comparison of carbon and nitrogen uptakes with 
(left) and without (right) the DOM model: (a) New, (b) 
regenerated (mol N m' 2 yr -1 ), (c) total carbon primary 
production (g C m' 2 yr’ 1 ) and (d) carbon excess uptake 
(g C m' 2 yr' 1 ) estimated by the DOM model for 2004 in the 
MAB and Georges Bank regions, (e), (f) and (g) are the same 
model outputs than (a), (b), and (c) for a simulation without 
the DOM model. Isobaths are in meters. 
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water column. The burial of PON has the same geographical distribution as POC with a C to N ratio 
of 9.3 (see Appendix). 


(a) station 1 
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(b) station 2 
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(c) station 3 
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Figure 8. Simulated daily primary production with (solid) and without the DOM (dots) and percentage of 
extracellular release (PER, dash) at (a) a 59 m depth station on Georges Bank (station 1), (b) a 55 m 
northern MAB station (station 2) and (c) a 46 m southern MAB station (station 3), 


4. Discussion 

4.1 Impacts of the DOM on the ecosystem model 

In the simulation without the DOM, the POM pools are directly remineralized to DIC and 
ammonium using the same rates as in Fennel et al. (2006), i.e. 0.03 d" 1 for small detritus and 0.01 d" 1 
for large detritus. Otherwise, the parameterization is the same for both simulations. The main impact 
of the introduction of the DOM is a large increase of regenerated production from 30% at Georges 
Bank and shelf break, 50% over the slope and deep ocean, to 250-300% in the mid- and inner-shelf 
of the MAB (Figures 7b and If). In contrast, the nitrate uptake shows globally the same distribution 
and level. The total carbon production shows an increase of 0 to 60 g C m'“ yr" in the open ocean 

2 i 

and the deeper Georges Bank (between 40 and 100 m, station 1), among which 0 to 20 g C m'" yr' is 

91 . 

linked to the carbon excess uptake. The increase in productivity ranges from 60 to 90 g C m' yr' in 

9 1 ... 

the outer-shelf and the major part of the slope and from 90 to 180 g C m'“ yr' in the inner- and mid- 
shelf among which 20 to 35gCm'~yr‘ is related to carbon excess uptake. This enhanced 
production is mainly (65 to 100%) caused by the progressive mineralization of the semi-labile DON 
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in surface waters in summer and autumn 
(Figure lie). The productivity supported 
by DON during summer and autumn could 
explain the usual underestimation of 
models which do not include DON (e.g. 
Fennel 2006). Only 0 to 35% of this 
increase is related to the ‘extra’ production 
of carbohydrates depending on the level of 
nutrient depletion. The shelf area that 
shows the minimum increase in primary 
production between the two simulations 
occurs in the simulation without DOM in 
the Georges Bank area between the 40 and 
100 m isobaths (Figure 8, station 1). In 
that area, the constant supply of nutrient 
by the tidal-induced mixing reduces the 
importance of DON as a source of 
nitrogen in the upper layer from 
mineralization. The phytoplankton 
biomass and chlorophyll levels are also 
dramatically lower without the DOM 
module, particularly in the MAB shelf and 
slope (Figure 6). This figure emphasizes 
accordingly a better estimate of the surface 
chlorophyll content in summer with the 
DOM components and higher levels in 
winter and autumn. 


4.2 The semi-labile DOM 
dynamics 

The most important contributor to the 
semi-labile DOC near the surface at 
station 3 is the phytoplankton exudation 
which occurs mainly within the mixed layer during the stationary and decaying phases of the bloom 
and between the mixed layer depth and the 10% isolume during summer (see Figure 5 and Figure 
12a). Note that the contribution of the carbon excess-based release in summer occurs in the upper 
part of this subsurface layer where phytoplankton is nutrient-limited (Figure 5c) and the nutrient- 
based release in the deeper and light-limited part (Figures 5a and 5b). The other important 
contributor to DOC release is the POC solubilization which occurs deeper in the water column than 
the phytoplankton release. The integrated flux of DOC release by POC solubilization (Figure 12c) is, 
however, about three times greater than the release by exudation although it is more diffuse within 
the water column. A large release occurs after the spring bloom below the 10% isolume due to the 
aggregation into fast-sinking particles. During summer, the release of semi-labile DOC by 
solubilization increases substantially between the surface mixed-layer and the depth of the 10% 
isolume where the small, slow-sinking, detritus dominates. The next increase of semi-labile DOC 
release by solubilization in summer occurs near the seabed in correlation with the production of 


Vertically integrated velocity (max =1.11 m.s -1 ) 



Figure 9. (a) Net horizontal transport of semi-labile DOC 
(mol C m' 2 yr -1 ) estimated by the model for the year 2004 in the 
MAB and Georges Bank regions: negative values are areas of 
production and export of semi-labile DOC and positive values are 
area of import, (b) Carbon burial (mol C m' 2 yr -1 ) for the same 
area and period. Isobaths are in meters. 
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Figure 10. Seasonal surface concentration of semi-labile DOC (mmol C m' 3 ) estimated by the model for 
the year 2004 in the MAB and Georges Bank regions. Isobaths are in meters. 



larger particles (fecal pellets) and the accumulation by POC resuspension. The smallest contributor 
to DOC release at station 3 is sloppy feeding by zooplankton (Figure 12b), which accounts for 
approximately 1% of the annual release, although it temporarily reaches 10% of the total DOC 
release. This simulated rate can reach 50% in highly productive areas such as the Chesapeake Bay 
mouth during a bloom in agreement with the field estimates of Moller et al. (2003). 

Overall at station 3, the vertically integrated reservoir of semi-labile DOC is 1.5 times higher than 
the carbon detritus pool and twice the carbon standing stock of phytoplankton (Figure 12d). At 
station 1 and 2, the semi-labile DOC pool is twice the POC pool and four times the carbon 
phytoplankton pool. It represents therefore the largest freshly produced organic pool of carbon in the 
water column of the mid and outer shelf. It can be efficiently exported by horizontal transport (see 
next section) since this organic carbon is in the dissolved form and is slowly mineralized. Figure 12d 
suggests that a large fraction of the carbon dioxide entering the shelf ocean is stored in semi-labile 
DOC since 12 to 14% of the C0 2 flux is buried at stations 2 and 3 (less at station 1 where tidal 
mixing prevents from deposition). 

4.3 Carbon export 

The annual carbon dioxide air-sea flux simulated by the model is positive for the entire study area 
(such as at station 3, Figure 12d) except in the upper part of the Chesapeake Bay where a continuous 


23 




12 Apr 16 May 19Jun 23 Jul 26 Aug 29 Sep 3 Nov 7 Dec 
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Figure 11. Seasonal profiles of (a) phytoplankton biomass (mmol N m' 3 ), (b) semi-labile DOC (mmol Cm' 3 ) 
and (c) semi-labile DON (mmol N m‘ 3 ) simulated by the model at station 3 (southern mid-shelf MAB). The 
mixed layer depth (solid black) and the depth of the 10% (dash white) and 1% (solid white) isolume are also 
presented. 


(a) Phytoplankton [mmolN m -3 ] at station 3 
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(b) Semi-labile DOC [mmolC m -3 ] at station 3 
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flow of terrestrial organic materials favors the mineralization and generates a source flux of 
dissolved CO 2 to the atmosphere. The annual flux is from 0.5 to 2.0 mol C m"~ yr" in the deep ocean 
and on the shelf except South of Cape Cod and Georges Bank where values range between 2 and 
5 mol C m’* - yr" . This is in agreement with DeGrandpre et al. (2002) who estimated a net annual 
uptake of ~1 mol m"“ yr" CO 2 on the MAB. Since the carbon-rich DOM buildup contributes to CO 2 
drawdown seasonally (Sambrotto et al., 1993), the results suggest that a large amount of the carbon 
entering the surface ocean is temporarily stored in DOC. 

Both degradation rates of the particulate and dissolved organic carbon (carbon detritus and DOC) are 
significantly lower than for the nitrogen detritus and DON leading to a lower regeneration rate for 
carbon than nitrogen. DOC and POC are consequently more efficiently exported through horizontal 
transport and sinking, respectively, than DON and PON. 

The carbon export of POC from the shelf to the slope has been studied extensively (e.g. Biscaye and 
Anderson, 1994; Thomas et al., 2002) and was shown to be particularly important near Cape 
Hatteras where both the MAB (Mayer et al., 2002) and SAB production (Schaff et al., 1992) 
contribute to the shelf-slope carbon efflux due to the converging shelf circulation. The comparison 
presented in Figure 9 shows that POC is buried in the inner- and mid-shelf of the MAB at rates 
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(a) Semi-labile DOC produced by exudation [mmolC m -3 d -1 ] at station 3 
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(b) Semi-labile DOC produced by sloppy feeding [mmolC m -3 d -1 ] at station 3 
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(c) Semi-labile DOC produced by solubilization [mmolC m" 3 d -1 ] at station 3 
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Figure 12. Seasonal profiles of semi-labile DOC fluxes released by (a) phytoplankton exudation 
(mmol N m' 3 ), (b) sloppy feeding and (c) POC solubilization (mmol C m' 3 d -1 ) simulated by the model at 
station 3 (southern mid-shelf MAB). (d) Vertically integrated pools of organic carbon and accumulated air- 
sea C02 and POC burial fluxes (mmol C m' 2 ) for the same station. The mixed layer depth (solid black) and 
the depth of the 10% (dash white) and 1% (solid white) isolume are also presented. 
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(d) Carbon pools vertically integrated or cumulated flux at station 3 (depth = 46 m) 
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comparable to the export of seasonally produced DOC from the outer-shelf and slope to the open 
ocean. In contrast to the southern MAB, the DOC produced at Georges Bank and south of Cape Cod 
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is mostly exported southward to the central MAB shelf between Long Island and Delaware Bay and 
does not contribute to a net export to the deep ocean. 

Conclusions 

The present study couples a circulation model to a biogeochemical model of carbon and nitrogen 
which includes the major DOM production processes to estimate (1) the role of DOM in the coastal 
ecosystem C and N cycling and (2) the relative importance of the export of freshly produced DOC to 
the open ocean compared to POC burial on the shelf. The carbohydrate production by 
phytoplankton, which occurs in nutrient-depleted and light-replete conditions, partially decouples the 
carbon and nitrogen primary production. The results show that DOM increases primary productivity 
by 60 to 180 g C m" yr' in the MAB, of which 65 to 100% is related to the ammonium release by 
DON mineralization in the upper layer and 0 to 35% to the ‘extra’ production of carbohydrates. 
Because DOM release by exudation and progressive mineralization occur near the surface, it is the 
most important process involved in this increase of primary productivity. However, in terms of flux, 
the annual release of semi-labile DOC by the deeper POC solubilization can be three times higher 
than the near-surface release by exudation. The seasonally produced DOC export from the shelf to 
the open ocean occurs mostly in the southern outer-shelf and slope of the MAB at a comparable rate 
(~1 to 2molm""yr") to POC burial in the inner- and mid-shelf. Subsequent steps in model 
development will consider the inclusion of the refractory DOC (as a passive tracer), multiple 
phytoplankton and zooplankton functional groups, a diagenetic sub-model to simulate 
remineralization and burial in the sediment, a fast sinking detritus (-100 m d" 1 ) and a higher 
horizontal resolution in shallow areas. The model parameterization and evaluation will also be 
improved by using new products derived from satellite remote sensing of surface DOC and POC 
concentration. These refinements will allow for a more complete estimate of the carbon budget at the 
scale of the Eastern U.S. continental shelf and provide for a better understanding of the role of DOC 
in the dynamics of carbon cycling at the land-ocean interface. 
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Appendix 


Parameter set 

See Tables Al, A2 arid A3 

Equations of the state variables 

Semi-labile DON and DOC 

The time rate of change of the semi-labile DON and DOC are: 

dDON/dt = Phytoplankton exudation + Sloppy feeding + 

+ Solubilization small and large N detritus (semi-labile fraction) 

- Remineralization semi-labile DON 

dDON 

— = M Phy + (1 - P) Qdon g Zo ° + 

dt 

- + ( s sDetN SDetN + s LDetN LDetN ) - a N 0 e ° 07 1 DON 

where p is the phytoplankton growth rate. 

dDOC/dt = Phytoplankton exudation (nutrient-based and carbon excess-based) 

+ Sloppy feeding + Solubilization small and large C detritus (semi-labile fraction) 

- Remineralization semi-labile DOC 


8DOC 

dt 


CN p (^ N // Phy + <r c y[ l i au L L (1 -L N )\Phy + ( 1 - /?) Q DOC S c g Zoo )+ 


- + SDetN SDetC + s LDetN LDetC ) - a co e 0 01T DOC ) 


where Ll and Ln are the non-dimensional terms that determine light- and nutrient-limitation, and 
Pmax the maximum phytoplankton growth rate (ju = ju max Ll Ln). 

Phy topi an kton 

Two sink terms are added in the phytoplankton time rate of change: the exudation terms of semi- 
labile and labile DON towards DON and ammonium respectively. 

dPhy/dt = Phytoplankton growth - Exudation of semi-labile DON 

- Exudation of labile DON (to NH4) - Grazing - Phytoplankton mortality 

- Aggregation with small N detritus - Sinking of living cells 

= // Phy (l — e n - co N ) - g Zoo - trip Phy - r ( SDet + Phv )Phv - w p 
dt dz 


where m p is the phytoplankton mortality rate, r the aggregation parameter of the small detritus and 
Phy (towards the large detritus pool) and Wp is the sinking velocity of living phytoplankton cells. 
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Phy is expressed in nitrogen unit using the constant C to N ratio ( CN P ) for accessing carbon units and 
therefore no equation is required for Phy expressed in carbon. A fraction ( &c ) of the carbon excess 
uptake represents the semi-labile DOC exudation by phytoplankton and is directed towards the semi- 
labile DOC. 


Chlorophyll 

The chlorophyll equation is modified accordingly to the changes of the phytoplankton equation: 

dChl/dt = Chlorophyll production - Loss by exudation of semi-labile DON 

- Loss by exudation of labile DON (to NH4) - Loss by grazing 

- Loss by phytoplankton mortality - Loss by aggregation with small N detritus 

- Loss by sinking of living cells 

— ^ = p chl ju Chl( 1 - s N - co N ) - g Zoo ^h — m p Ch] - r(SDetN + Phv)ChI - w p L^JlL 
dt 1 1 Phy dz 


where p c hi is the fraction of phytoplankton growth devoted to the chlorophyll synthesis (Geider et al., 
1997): 

Pchi = o max p Phy / al Chi 

where 6 max is the maximum ratio of chlorophyll to phytoplankton biomass, a is the initial slope of 
the phytoplankton growth curve relative to light and I the photosynthetically available radiation. 

Zooplankton 

The zooplankton, like the phytoplankton, is only expressed in nitrogen unit: 

dZoo/dX = Fraction of grazing assimilated 

- Excretion (basal metabolism and grazing dependent) - Mortality 


dZoo 

dt 


= g fiZoo - 


+i E /?-§- 

g max 


Zoo-m z Zoo~ 


T\vq remaining term [(l-f3) g Zoo] is divided between the production of semi-labile and labile DON 
by sloppy feeding (towards semi-labile DON and DIC respectively) and the production of fecal 
pellets (towards the small N detritus pool). 

The zooplankton equation expressed in carbon is: 


dZoo 

dt 


I c 


CN p g p Zoo - CN Z 


f 

^BM 

V 


+ 1 E /? 


g 

gmax 


Zoo - m z Zoo~ 

/ 


-r Cac CH z Zoo 


where l'cexc is the rate of carbon excess respiration due to the C to N ratio difference between 
phytoplankton and zooplankton. The constant zooplankton C to N ratio ( CNz =5.0) leads to the 
formulation: 


CN Z = (CNp p g Zoo - r Cexc CN Z Zoo)/ p g Zoo 
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Or 


I’C'exc pg(CN P -CNz)/CN z 

This excess of respired organic carbon is directed to DIC. It ensures the conservation of the 
zooplankton C to N ratio and therefore the zooplankton equation expressed in carbon is implicit. 

DIC 

The air-sea exchange of carbon dioxide is taken from Fennel et al. (in prep.). 

dDIC/dt = - Nutrient-based uptake by phytoplankton growth 
- C excess-based semi-labile DOC exudation 
+ Nutrient-based exudation of labile DOC 
+ Labile DOC produced by sloppy feeding 
+ Excretion (basal metabolism and grazing dependent) 

+ Solubilization small and large detritus C (labile fraction) 

+ Remineralization of semi-labile DOC 
+ Air-sea C02 flux 


dDIC 


dt 


= CN P (- fuPhy - ycr c /J max L L (1 - L N )Phy + co c /uPhv + (1 - (3)Q DOC ( 1 -S c )gZoo ) + 


( 


... + CN „ 


■■■+-— f^-C0 2 (pCO atr - pC()) 
bz 


g 




^ r Cexc 

g max J 


Zoo + S SDelc SDetC + S LDttc LDetC + a co e° mT DOC + 


Ammonium 

dNH 4 /dt = - Ammonium uptake by phytoplankton growth + Exudation of labile DON 
+ Labile DON produced by sloppy feeding 
+ Excretion (basal metabolism and grazing dependent) 

+ Solubilization small and large detritus N (labile fraction) 

+ Remineralization of semi-labile DON - Nitrification 


dNH4 ( g 

— — = -M ma L L L NH4 Phy+co N fjPhy + (l-fi)Q DON (l-S N )gZoo+ l BM +l E p— 
dt l g 


\Zoo + 


V 


max / 


... + (1 - 5 N ){S SDem SDetN + S w m LDetN) + a m e wlT DON - nNH4 


where Ibm and Ip are the zooplankton excretion rates due to basal metabolism and assimilation 
intensity respectively, and n is the nitrification rate (same parameterization than in Fennel et al. 
(2006). Ll is the non-dimensional light limitation and Lnh 4 is th e nutrient limitation term for 
ammonium. 
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Nitrate 


dN0 3 /dt = - Nitrate uptake by phytoplankton growth + Nitrification 
'' W! = ~M^L L L Nm Phy + n NH4 


where Lno 3 is the nutrient limitation term for nitrate. 

Detritus N 

dSDetN/dt = Phytoplankton mortality - Aggregation with living phytoplankton cells 
- Small detritus N solubilization - Sinking of small detritus N 


DSDelN = m ^ phy _ T ( SDetN + phy'jSDetN - s SDelN SDelN - w> s 5SDetN 


dt 


dz 


dLDetN/dt = Fecal pellets production (nitrogen fraction) + Zooplankton mortality 
+ Aggregation of small detritus (N) and phytoplankton cells 
- Large detritus N solubilization - Sinking of large detritus (N) 


= (1 - P)(\ - Q DON )gZoo + m : Zoo 2 + t (SDetN+ Phy f - W LDetN - w L dLDelN 
ot dz 

Detritus C 

5SDetC/^t = Phytoplankton mortality (C) - Aggregation with living phytoplankton cells 
- Small detritus C remineralization - Sinking of small detritus C 

5SDerC = CN p m p Phy - z (SDetC + CN p Phv)SDetC - s SDelc SDetC - 

dt dz 

dPOCi/dt = Fecal pellets production (carbon) + Zooplankton mortality (C) 

+ Aggregation of small detritus (C) and phytoplankton cells (C) 

- Large detritus C solubilization - Sinking of large detritus (C) 


dLDetC 

dt 


= CN P (1 - p)(\ - Q doc )gZoo + CN z m ; Zoo 2 + z [SDetC + CN p Phyf - s LMC LDetC + 


... - w 


dLDetC 
L dz 


where ws and wi are the sinking velocities of small and large detritus respectively. 


Oxygen 

The oxygen equation taken from Fennel et al. (in prep.) is modified as the following: 


30 



dOx/di = Phytoplankton growth linked to N uptake 

+ Phytoplankton growth linked to C excess uptake 

- Labile DOM oxidation - Labile DOM oxidation from sloppy feeding 

- Zooplankton excretion 

+ Solubilization small and large POM (labile fraction) - DOM oxidation 

- Nitrification +Air-sea flux 


dOx 

dt 


(Av03 r 02:N03 F NHA ^Ol.NHA )/ / ma x F L F ^y & c CN p /^ max ^ z (1 L N )PflV + 


... - K 


02.NH A 


co N Phy + (1 - P)Q don {\ - S c )gZoo + 

1 BM 1 eP~ 


g 




...+ (1 - S N )(s SDetN SDetN + s LDetN LDetN) + a m e 0 07T DON_ 


S max ) 


Zoo + 


+ 


vK 

... - 2 n NH 4 + — —{Ox sal - Ox) 
A z 


where r 0 2 :N 03 and r 0 2 :NH 4 are the 0 2 :N ratio for nitrate and ammonium respectively, vK 0 2 is the gas 
exchange coefficient for oxygen, Az is the thickness of the top box and Ox sa t is the saturation 
concentration of oxygen. 

The oxygen produced by the synthesis of carbohydrates (carbon excess uptake) has a one to one 
mole ratio with DIC following the equation: CO 2 +LL0 + energy -> (CLfO) + O 2 . 

The oxygen uptake from the solubilization of POM (PON+POC) and the oxidation of DOM 
(DON+DOC) is approximated to the solubilization of detritus N and oxidation of DON with a C to 
N ratio of 6.6. 


Bottom boundary condition 

In order to take into account the resuspension of detritus C near the seabed due to bottom friction, a 
fraction (A, ies , see next section) of the bottom carbon flux is resuspended and mineralized in the lower 
water column. The complementary fraction (1-kres) is buried, the flux of buried carbon thus is: 

F = BE (\-l )F 

1 Cbwied '\es) L Cbottom 

Where BEc is the burial efficiency (see next section) and F cbottom is the detritus C flux that reaches 
the bottom before resuspension: 


F =CN w 

1 Cbottom ^ y P VV P 


dPhy 

dz 


+w c 


I z=H 


dSDetC 

dz 


+w r 


:=H 


dLDetC 
dz 


I z=H 


The POM which is not resuspended nor buried is mineralized and therefore the bottom boundary 
condition for carbon follows: 

dDIC/dt| z =H = mineralization of (resuspended bottom detritus C 

+ not resuspended nor buried bottom POC) 


dDIC 

dt \ 


~ F Cbottom fares + ( l ~ Kk )(1 ~ BE C )) 


z=H 
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For nitrogen, the same resuspension rate is applied (A res ) to the detritus N reaching the seabed. The 
remaining detritus N is subject to burial and denitrification following: 


F = ; F 

Nresuspended res Nbottom 
E Nburied ~ BE n (1 — X res ) F } 


Nbottom 


R 


Ndeni trifled 


= (\-BE N )(\-X res )F, 


Nbottom 


where F Nbottom is the flux of detritus in nitrogen that reaches the bottom before resuspension: 


F, 


Nbottom 


= W, 


dPhy 

dz 


+w c 


dSDetN 

dz 


+ W T 


dLDetN 

dz 


The denitrification process is taken into account as it was shown to be significant in the MAB 
(Fennel et al., 2006). The stoichiometry calculation shows that the bottom boundary condition for 
ammonium is: 


dNFLf dt|z=H = mineralization of (resuspended bottom detritus N 
+ not resuspended nor buried bottom detritus N 
+ not resuspended nor denitrified bottom detritus N) 


dNH A 

dt 

dNH A 

dt 


I z=H 


= F + — p 

Nresuspended , r Ndenibified 

to 


Z=H V 


'L+t- O-'U (\~be n ) 

f 6 


F 


Nbottom 


The total amount of nitrogen lost through burial and denitrification (N 2 ) is: 


F - F +— F 

Nlost Nburied , r Ndemtnfied 

to 


E yiost ~ , r ^res + EE N )F mottom 

f 6 


The POM resuspension 

The resuspension is taken into account as a function of the friction velocity at the seabed (U ). The 
resuspension rate (%) follows: 



^ , m m • m • • • •••• 

where U d is the critical friction velocity above which all organic matter is maintained in suspension 
(U </ = 0.31 cm.s’ , Peterson, 1999). The resuspended fraction of POC is thus largely dependent of 
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the local near bottom current velocity which is driven by the general circulation and the tides on the 
continental shelf and also by wind events in shallow waters. 

The POM burial 

Thomas et al. (2002) measured and reviewed high rates of carbon burial along the continental shelf 
of the U.S. northeastern continental shelf. The carbon and nitrogen burial rates have been 
implemented to simulate the loss of material in the sediment. 

A fraction (burial efficiency, BEc) of the particulate organic carbon that reaches the seabed is buried 
following the empirical expression of Henrichs and Reeburgh, 1987: 

log Fc = 0.69 log w + 2.27 

and 


BE C = w 4 72.1 

where Fc is the organic carbon flux at the sediment surface (gC m" y" ) and w is the sediment 
accumulation rate (cm y" 1 ). Resolving the system leads to the following formulation for BE C (%): 


BE 


c 


1 

2.1 


10 


loe Fc\ 

— -2.27 

0.69 ) 


-i 0.4 


This formulation matches the upper values of Thomas et al. (2002) who measured and reviewed 
burial efficiency values (in % of organic carbon deposition): 10-20% at the slope off Cape Cod 
(SEEP-I), 25-50% in the MAB (SEEP-II) = 25-50, 3-40% at the slope off Cape Hatteras. 

Gelinas et al. (2001) reviewed the C to N ratio of buried organic matter and reported values of 9-10 
for the shelf and estuarine surface sediments and slightly lower in deeper waters. A value of 
CNbunai= 9.3 is used to estimate the flux of buried organic nitrogen in agreement with values 
measured in the sediment of the MAB shelf (Mayer et al., 2002; Mayer, unpub. data). 

A maximum of 75% of carbon burial efficiency is applied as it corresponds to the maximum value 
measured: 

BE C = MIN{ [1 0(log Fc 70.69-2.2 7) ]0. 4/2. 1 ; 0. 75 } 

For nitrogen burial, a similar expression of burial efficiency is used introducing a CNburiai ratio: 

BE n = MIN{ [10 (I ° 8( CNburialFN)/ 069 ~ 2J7) ] 04 72.1 ; 0.75 } 


33 



Tables 


Table 1: Summary of the terms involved in the uptake and exudation of nitrogen and carbon by 
phytoplankton. PPl is the nutrient-based primary production limited by light, Ln is the nutrient 
limitation, CNp is the C to N ratio for phytoplankton and y the parameter of carbon excess-based 
DOC excretion by phytoplankton (see Table A1 for the definition of other parameters). 


Expression 


Description 


Ui v = PPl Tjv 
Unut.C = CNp U N 

E lab.N = C0 N Un 
E lab. nut. C = CNp CO N U N 
E sem. N = £iV U N 
Esem.nut.C CNp E sem jy 

Uexc.c = CNp y PPl (1— L n ) 

Ejab.exc.C (1 @c) C e:xc c 

Esem.exc.C @c Uexc.C 


nitrogen-based primary production or uptake of nitrogen 
nutrient-based primary production in carbon 
exudation of labile DON (directed to ammonium) 
exudation of labile DOC (directed to DIC) 
exudation of semi-labile DON 
nutrient-based exudation of semi-labile DOC 
carbon excess uptake by nutrient-stressed phytoplankton 
carbon excess-based exudation of labile DOC 
carbon excess-based exudation of semi-labile DOC 


Table 2: Release DCAA flux from aggregates, carbon content of aggregates, fraction of nitrogen in 
released DCAA estimated by Smith et al., 1992), and DON release rates deduced from thes e data. 


Aggregate 
type (April 
1990) 

DCAA release 
(pg.aggfd' 1 ) 

(1) 

Carbon content 
(pg C.agg" 1 ) 

(2) 

Nitrogen fraction in 
released DCAA 

(non-dimensional) 

(3) 

DON release rate 

( DC AA • §N ag g .d ) 

assuming POC:PON=5 
(l)*5.0*(3)/(2) 

Larvacean 

house 

0.936 

3.5 

23.3/155.2=0.150 

0.200 

Diatom floe 

0.527 

3.5 

2.4/15.4=0.156 

0.117 

Larvacean 

house 

0.478 

3.2 

8.0/53.6=0.149 

0.118 

Larvacean 

house 

0.365 

4.5 

0.5/3.1=0.161 

0.066 
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Table Al: DOM specific parameters 

Symbol Value or range and Unit 

co N 0.03 [nondimensional] of 

primary production (N) 
co c 0.03 [nondimensional] of 

primary production (C) 
s N 0.04 [nondimensional] of N 

primary production 
y 0.20 [nondimensional] 

o c 0.45 [nondimensional] 

S N 0.30 [nondimensional] 

S c 0.55 [nondimensional] 


0.0-0.71 [nondimensional]: 

function of the ratio 

grazing: maximum grazing 

(§/ §max) 

0.0-0.71 [nondimensional]: 

function of the ratio 


grazing: maximum 

(§/§max) 

0.01530 d' 1 


grazing 


a co 

0.00767 d' 

s SDetN 

0.11 d’ 1 

$LDetN 

0.11 d’ 1 

s SDetC 

0.08 d’ 1 

SLDetC 

0.08 d’ 1 


Parameter or formulation 

Labile DON exudation rate 

Nutrient-based labile DOC exudation rate 

Exudation rate of phytoplankton semi-labile DON 

Parameter of carbon excess-based DOC exudation 
Fraction of semi-labile DOC produced by the carbon 
excess-based exudation 

Fraction of semi-labile DON to total DON within the 
phytoplankton cell 

Fraction of semi-labile DOC to total DOC within the 
phytoplankton cell 

Fraction of total DON to (DON+PON) within the 
phytoplankton cell 

Fraction of total DOC to (DOC+POC) within the 
phytoplankton cell 


Remineralization rate of semi-labile DON at 0°C (a NT 
= a N0 e°' 07T , with T in °C) 

Remineralization rate of semi-labile DOC at 0°C (a CT 

= a co e 0 07T , with T in °C) 

Bacterial solubilization rate of small N detritus 
Bacterial solubilization rate of large N detritus 
Bacterial solubilization rate of small C detritus 
Bacterial solubilization rate of large C detritus 
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Table A2: Modified parameterization from Fennel et al. (2006) 


Symbol 

New value or range, 
and Unit 

Former value 

Parameter 

a a 

0.020 (W.m-2 ) _1 . d' 1 

0.025 (W.m-2 )'\ d' 1 

Initial slope of the P-I curve 

a 

Mmax 

1.6 d 1 

0.59*1.066 T (J is the 
temperature in °C, 
Eppley, 1972) 

Maximum growth rate of 
phytoplankton 

CN Z b 

5.0 [nondimensional] 

6.625 

Zooplankton C to N ratio 

P b 

0.22-085 for 

g/gmax=[ E 0-0.0] 

[nondimensional] 

0.75 

Zooplankton assimilation 

efficiency 

c 

w s 

1.0 m d’ 1 

0.1 md' 1 

Small detritus sinking velocity 

c 

W L 

10.0 md' 1 

1.0 m d' 1 

Large detritus sinking velocity 


a Brush et al. (2002) showed that the temperature-dependent formulation of Eppley (1972) 
underestimates primary production. Even if a temperature dependency most probably exists in 
relation to the cell metabolism, the light intensity is the prior control factor. The use of a 
temperature-dependent formulation led to a latitudinal variation and underestimation (low 
temperature below the thermocline) of primary production. For a better analysis of the results, the 
temperature dependency of the maximum growth rate was totally removed. 
b See text for details. 

c Since the dead phytoplankton cells on one hand, and the zooplankton corps and fecal pellets on 
the other hand sink with distinct velocities due to their particles size difference, the zooplankton 
products are flowed to the large particle pool which sinks faster instead of the small particle pool. 
The aggregation process thus concerns only the phytoplankton living and dead cells. The sinking 
velocities proposed for such a configuration are 1 m d' 1 for the small detritus pool (dead 
phytoplankton cells) and lOmd' 1 for the large detritus pool (zooplankton particles and 
phytoplankton aggregates). 


Table A3: Common parameterization with Fennel et al. (2006) 


Symbol 

Value and Unit 

Parameter 

^NOS 

0.5 mmol N m' 3 

Half-saturation constant for nitrate uptake 

k NH4 

0.5 mmol N m’ 3 

Half-saturation constant for ammonium uptake 

CNp 

6.625 [nondimensional] 

Phytoplankton C to N ratio 

Smax 

0.6 d' 1 

Maximum grazing rate 

Kphy 

2.0 (mmol N m' 3 ) 2 

Half-saturation constant for grazing 

Wp 

0.15 d' 1 

Phytoplankton mortality 

T 

0.005 (mmol N m' 3 )' 1 d' 1 

Aggregation parameter 

Qnax 

0.053 mgChlfl mgC" 1 

Maximum chlorophyll to phytoplankton ratio 

Ibm 

0.1 d' 1 

Excretion rate due to basal metabolism 

Ie 

0.1 d' 1 

Maximum rate of assimilation-related excretion 

m z 

Wphv 

0.025 (mmol N nf 3 )' 1 d' 1 
0.1 m d' 1 

Zooplankton mortality 
Phytoplankton sinking velocity 
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